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p^ . The goal of combining beamforming and space-time coding in this work is to obtain full-diversity 

Cn . order and to provide additional received power (array gain) compared to conventional space-time codes. 

In our system, we consider a quasi-static fading environment and we incorporate both high-rate and low- 
rate feedback channels with possible feedback errors. To utilize feedback information, a class of code 
c/3 i constellations is proposed, inspired from orthogonal designs and preceded space-time block codes, which 

IS called generalized partly orthogonal designs or generalized PODs. Furthermore, to model feedback 
errors, we assume that the feedback bits go through binary symmetric channels (BSCs). Two cases are 
^ I studied: first, when the BSC bit error probability is known a priori to the transmission ends and second, 

^.f-N I when it is not known exactly. In the first case, we derive a minimum pairwise error probabiHty (PEP) 

'^ ' design criterion for generalized PODs. Then we design the quantizer for the erroneous feedback channel 

• ' and the precoder codebook of PODs based on this criterion. The quantization scheme in our system 

f^ ' IS a channel optimized vector quantizer (COVQ). In the second case, the design of the quantizer and 

the precoder codebook is based on similar approaches, however with a worst-case design strategy. The 
attractive property of our combining scheme is that it converges to conventional space-time coding with 
low-rate and erroneous feedback and to directional beamforming with high-rate and error-free feedback. 
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}_j ■ This scheme shows desirable robustness against feedback channel modeling mismatch 
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I. Introduction 

We study point-to-point closed-loop communications over multiple-input single-output (MISO) 
quasi-static Rayleigh fading channels. This channel model characterizes the downlink of several 
wireless communication systems such as the new generation of fixed and mobile cellular systems 
in IEEE 802.16 [1]. We assume that the single-antenna receiver is able to estimate the channel 
coefficients perfectly. Acquiring perfect receiver channel state information (CSI) is relatively 
easy. Because by definition, in quasi-static channel environments, there is a long fading block 
length and accurate channel estimation through training is possible. Obtaining transmitter CSI 
(CSIT), however, requires the use of feedback in the system. When CSIT is available either 
fully or partially, directional beamforming can increase the average received signal-to-noise ratio 
(SNR). This benefit of beamforming is called array gain. It improves the capacity of wireless 
communication systems [2]-[4], reduces the outage probability [5] [6], and enhances the error 
performance. Partial CSIT is widely used in the literature for combining space-time coding and 
beamforming using precoded space-time codes, [1], [7]- [11], and for combining beamforming 
and power control [12] [13]. 

We design a fixed-rate practical coding and beamforming system with instantaneous power 
constraint at the transmitter. The performance measure is the pairwise error probability (PEP) 
of transmission over a finite length data frame. The diversity order in this context is defined as 
the decay rate of the pairwise error probability with SNR. In MISO systems, full-diversity order 
means an error probability decay rate equal to the number of transmit antennas. In quasi-static 
environments, channel fading coefficients remain constant during the transmission of each data 
frame, which occurs within the channel coherence time, and change independently afterwards. 
Therefore, feedback information must be updated in every frame. Usually, feedback channels 
are severely bandwidth and power limited. It is preferable to reduce the feedback rate, using 
resolution-constrained (quantized) CSI [14]- [19]. Most of the work in the literature assume error- 
free (noiseless) feedback. However, this assumption requires tremendous protection of feedback 
information in a practical system, which requires the dedication of excessive bandwidth and power 
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to feedback channels. We study a practically motivated scenario, where not only is the feedback 
link resolution constrained, but also the existence of noise in the feedback link introduces errors 
to feedback information [20]. This situation can naturally arise in a variety of applications, where 
feedback information must be carried through actual fading channels [21]. If the system is designed 
robustly against feedback errors, a significant portion of the additional bandwidth and excessive 
power spent for protecting feedback information can be saved. This can increase the efficiency of 
the overall network. 

For simplicity, we consider a simple discrete memoryless channel (DMC) model of the feedback 
link, using binary symmetric channel (BSC) models for feedback bits. We assume that the cross- 
over probability (bit error probability) of feedback bits is fixed over all channel realizations. This is 
a simple, yet practical version of the bit error model in [22], which introduces a finite- state fading 
channel model. A more elaborate model of feedback errors can be found in [23]. Other than the 
limited bandwidth and feedback errors, possible delay in feedback can also be problematic. With 
delay, feedback information might become outdated. Therefore, the system performance degrades 
due to the mismatch between the CSIT and the real channel realization [13] [24]. In this work, 
however, we only consider a simple erroneous feedback channel model to establish the design 
guidelines and leave the extensions to other feedback channel models for future work. 

If feedback information is error-free, beamforming or combining space-time coding and beam- 
forming can provide full-diversity order on top of array gain [25] [26]. However, with feedback 
errors, achieving full-diversity order is not straightforward [20] [21]. It was recently shown that 
with less restriction on the feedback rate, we can obtain full-diversity order by combining space- 
time coding and beamforming [25]- [32]. This is even possible with feedback errors in some 
scenarios. For example in [33] this issue is addressed by appropriate spatial power allocation 
using mean and co variance feedback in block-fading channels. In our combining scheme in this 
work, the goal is to address the low-rate feedback and feedback error problems at the same time 
for quasi- static channels. Moreover, we want our scheme to extend easily for any feedback rate 
and feedback error probability, even if there is a mismatch between the actual feedback channel 
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parameters and the knowledge of the transmitter/receiver sides about these parameters. 

The rest of the paper is organized as follows: In Section II, the system model is sketched and 
our combining scheme is introduced. This scheme is based on a code structure named generalized 
partly orthogonal designs (PODs). Then we derive the design criteria of our combining scheme, 
based on PEP analysis. In Section III, using channel optimized vector quantizers (COVQs) similar 
to [35]- [38], we propose a quantizer structure for the feedback channel and a precoder structure 
for generalized PODs. First we consider the case when feedback error probability is known. The 
case of feedback channel modeling mismatch, or not knowing feedback error probability a priori 
is also discussed. The latter discussion is inspired from the worst-case COVQ design approaches 
in [39]. Section IV provides the numerical results of the paper and finally Section V draws our 
major conclusions. 

Notations: In the sequel, p(i), p(j|i), and p(j;i) denote the marginal probability of the event 
i, the conditional probability of the event j given i, and the joint probability of events j and 
i, respectively. The operators || • \\f, {-Y , {■)\ and (■)* represent Frobenius norm, transpose, 
Hermitian, and conjugate, respectively. f(x) denotes the probability density function (pdf) of the 
random vector x and f(x|z) denotes the pdf of x given event i. Finally, Ex{^(x)} shows the 
expectation of a function of x and Ex|j{\E'(x)} shows the same expectation given event i. 

II. System Components and Design Criteria 

The block diagram of our system is shown in Fig. [Tl In this system, the receiver is equipped 
with perfect channel estimation, ideal synchronization, and maximum-likelihood (ML) decoding. 

A. Forward and Feedback Channel Parameters 

During the transmission of each data frame, we represent the MISO quasi-static fading forward 
channel coefficients by an M-dimensional complex Gaussian vector h. The task of the quantizer 
is to quantize this vector at the receiver and generate a quantization index. The feedback channel 
conveys the quantization index back to the transmitter. Due to the limited feedback bandwidth, in 



some cases we quantize a portion of h by decomposing this vector into two parts as 

"(M-Af)xl ^%xl 



(1) 



where h is the A^-dimensional quantized part and h is the (M — A^) -dimensional part which 
remains unknown to the transmitter. The parameter N in our system that is the dimension of 
the quantized vector depends on the number of feedback indices in the feedback codebook. Note 
that we deal with short-term (instantaneous) power constraint at the transmitter. Also our scheme 
only supports fixed-rate transmissions. Therefore, the amplitude of the quantized channel vector is 
irrelevant to the performance and the design of our system. Our quantization scheme works on the 
direction of the channel sub-vector h, i.e., h = h/||h||i7. We decompose h into its amplitude and its 
direction through h = ^/7 h, where 7 is the amplitude square. The direction vector h is uniformly 
distributed over the unit amplitude complex sphere C^ and the amplitude square 7 is a Chi-square 
random variable with 2N degrees of freedom. The pdf of 7 is £(7) = 7^^^ exp (— 7)/(A^ ~ 1)' • 

In our feedback scheme, we partition the space of h into K mutually exclusive feedback 
(Voronoi) regions V = {Vi, ■ ■ ■ ,Vk}, whose union spans the whole unit-norm complex sphere 
C^. The objective of the quantizer at the receiver is to encode h by finding the index of the Voronoi 
region Vj that h belongs to. The objective at the transmitter is to choose a precoder matrix Vj 
based on the transmitter's feedback index j, which may be different from i as a result of errors 
in the feedback channel. We first need to design a precoder matrix codebook V = {Vi, ■ ■ ■ , Vk} 
and store it at both ends of the communication link. The precoder matrix Vj defines a matrix 
constellation set at the transmitter. 

In our system, after finding a feedback index, the receiver converts it to bits and sends the 
bits to the transmitter. The typical value of the feedback bit error probability is approximately 
0.04 [20] [21]. However, for the sake of generality, we will consider a case where feedback bit 
error, pj, can take any value between 0.00 and 0.50. Note that if pj > 0.50, the transmitter can 
simply flip the feedback bits and achieve better performance results. In the sequel we assume that 
the transmitter/receiver sides of our system know pf a priori, unless stated otherwise. Suppose 
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that the feedback index i E {1, ■ ■ ■ , K} is chosen at the receiver. The probability of this event 
is denoted by p(i). Also, suppose that the mapping of feedback indices to feedback bits is an 
identity mapping. For any feedback index i, the feedback bits can be simply obtained as binary 
representation of i — I. The encoder of the quantizer (and the mapper) are located at the receiver 
and the decoder of the quantizer (and the demapper) are implemented using a simple table look- 
up operation at the transmitter. In our system model, feedback indices go through a DMC from 
log2 K uses of BSCs with cross-over error probabilities pf for each feedback bit. The feedback 
channel converts i to j E {I, ■ ■ ■ , K}, with the conditional index inversion probability 

where d{i,j) is the Hamming distance between the binary representations of j and i. 

B. Modulation Scheme 

The transmitter uses a class of coded modulation schemes called generalized partly orthogonal 
designs (PODs). PODs were firstly introduced in [40]. Using PODs, we can utilize feedback 
information from any number of channel coefficients A^ for combining coding and beamforming 
across M transmit antennas, where M > N. The precoder or beamformer structures incorporated 
in original PODs of [40] were designed based on maximizing the received SNR. The generalized 
PODs of this work are however designed based on minimizing the pairwise error probability 
and hence they have a different structure. Unlike the original PODs, the precoder structures used 
in generalized PODs allow spatial power allocation, depending on the probability of error in the 
feedback link. Like original PODs, generalized PODs also have two parts: i) the STBC inner code 
that uses an M x T orthogonal design from [41] [42] and ii) the precoding part, which uses an 
N X N (A^-dimensional) precoder matrix, chosen based on feedback information. Each modulation 
matrix in the constellation set is limited to a short-term (instantaneous) power constraint MT. 

For clarification, let us use the following example. Suppose that the MISO forward channel 
has 4 transmit antennas. Similar to [40], we construct the code structures on the columns of an 
orthogonal design matrix. With a slight abuse of the notation, in this example we use the following 



M X T orthogonal design structure where M = T = 4 



Zl — Z2 — Z3 —2:4 

Z2 21 Z4 -Z3 

Z3 —Z4 Zi Z2 

Z4 Z3 -Z2 Zl 

Here, 2;^ : /t G {1, ■ ■ ■ ,T} denote the data symbols chosen from a real symbol constellation. 
Generalized PODs allow using any precoder dimension N, where A^ < M. The dimension of 
the precoder depends on the dimension of the channel vector quantizer. For instance, using a 
2-dimensional precoder, combining is performed through the following code structure: 
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Also, using a 4-dimensional precoder, we can use the following code: 
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In these codes, p^^^ denotes an A^xA^ precoder matrix with power N, i.e., Tr ( Vf^'^V, ^ 
N, where Tr(-) denotes the trace operation. Note that when N = M generalized PODs look similar 
to a conventional PSTBC in [25]. Also when M ^ N, they look like a POD in [40]. However, 
they are different structures, as we will clarify further in the sequel. Using similar approaches, 
one can design generalized PODs from complex orthogonal designs for any dimensions [42]. 
Furthermore, similar ideas can be applied to design a generalized partly quasi-orthogonal design 
based on quasi-orthogonal inner STBCs from [43] [44]. The main challenge is how to design the 
precoder matrix p^^^ to obtain minimum pairwise error probability with each code structure. 

Suppose that the quantizer codebook has K indices. We will show that in order to obtain full- 
diversity order using the above codes, we must choose precoder structures with dimension N < K. 
li K > M, we have high-rate feedback and PODs with N = M-dimensional precoders can provide 
full-diversity. Furthermore, we will show that in an M-dimensional MISO system, where K < M, 



which resembles low-rate feedback in our system, the least pairwise error probability is associated 
to the code structures with A^ = fT-dimensional precoders. As a result, to minimize PEP using a 
generalized POD, we should use N = min{i^, M} -dimensional precoders. 

C. Pairwise Error Probability Analysis 

Using generalized PODs, the equivalent base-band signal at the receiver can be modelled as 

y = Z/ h + n (5) 

where Z is the transmit signal matrix (as the Hermitian of a generalized POD modulation matrix) 
and n is the T-dimensional noise vector. Note that for simplicity we dropped the upper index of 
Zj. If the regulated average SNR at the receiver is t^q, each element of n is a complex circularly 
symmetric additive white Gaussian noise (AWGN) variable with variance c^ = l/(M?7o). 

To decode the matrix Zj, we assume that the index of the constellation set, j, which rep- 
resents the modulation scheme at the transmitter is known at the receiver. Usually, in wireless 
communications standards, the transmitter sends some control signals to the users (receivers) in 
the header field of each data frame that indicate the modulation scheme. For example in IEEE 
802.16, downlink burst profile of the physical layer, which is a part of downlink channel descriptor 
of the MAC layer contains the type of modulation used [1]. Our system requires to include the 
index of the constellation set, j, in this category and to assume that the receiver knowledge of j is 
updated as frequently as feedback is applied. By this assumption, the conditional PEP of detection 
in favor of an erroneous codeword Zj when Zj is transmitted can be tightly upper bounded by 
the following Chernoff bound [42] 

p(Z,^Zj|h)<-exp^-^^j (6) 

where D(Zj,Zj) = h^ (Zj — Z'j){Zj — Z'jf h. Now, suppose that the receiver processes fi and 
chooses Vj from the set of Voronoi regions V. The average probability of pairwise error given the 



feedback index i at the receiver can be expressed as 

p{z, ^z'^\i)= I ^ p(z, ^ z;|h) f(h| i) dh (7) 

^Dom(h|i) 

where Dom(h|z) is the domain of the random variable h, conditioned on the receiver's feedback 
index. The variable h can be expressed as a one-to-one function of the three variables, h, h, and 7. 
Therefore, f(h| i)d\i is statistically equivalent to f(h; 7; h\ i) dh d'j dh [34]. Also, for independent 
and identically distributed (i.i.d.) channel realizations, the latter three variables are independent. 
Therefore, f(h; h]'y\ i) = f(h) £(7) i(h\ i) and we can reexpress (|7]) as 

p(Zj- ^ z'j\i)= I [ [ p{Zj -^ z'j I h, h, 7) f(h) i(h I i) f(7) dh d-f dh (8) 

JDom(h) JtI+ JVi 

Here, 7^+ = [0, 00) is the domain of 7 and Dom(h) is the domain of the un-quantized portion of 
the channel vector. For generalized POD structures, using Equation Q, the above PEP expression 
can be upper bounded by: 



p{z^^zr\i) <]- f dh [ dj [ dhiih)ii^)m^) (^^vi-VcliP + ^]) (9) 

^ JDom(h) Jn+ Jv^ 

where (3 = h VjVjh and ric = \J2k.=i I-^k ~ ^'lifl/'^'^l- Note that rjc is proportional to the 
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takes values on 

F 



Euclidian distance of the inner code matrices and the SNR. Also ^ = 
[0, 00) and follows a Chi-square distribution. In the sequel, with a slight abuse of the notation, 
we denote pj^^^ by Vj. Because, the precoder dimension N is fixed and the transmitter must 
pick the precoder matrix from V = {Vi, ■ ■ ■ ,'Pk} using feedback index j. To proceed, we use 
the following relations: 



f(^) exp i-T],^) d^ = {l + Vc)'^^^'^^ 


exp (-r/c 7 P ) f(7)c?7 = (1 + VcP)'^ 
Then the PEP of the worst-case error event, conditioned on the receiver index i is 

p(Z,. -. Z] I < ^ (1 + r^c)"^^'"^^ %,v, { (1 + ^c^V.PJ^) '"" } (10) 
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where r]c is the minimum Euclidian distance between the inner STBC parts of Zj and Z^ divided 
by the noise power. Note that 1]^ represents the transmission signal power divided by the noise 
power at the receiver antenna or the transmit SNR. When we use unit-norm symbols such as 
BPSK, in a worst-case error event we have r]c = M?7o/(4T). 

D. Expanding PEP and Deriving the Design Criterion 

Our goal is to minimize the average PEP over all feedback realizations: 

K K 

i=i j=i 

K K 

= E Ep(^'' ^)pi2,^Z'^\J■^) (12) 

i=i j=i 

K K 

= E Y.'PiiM Pi^ Pi2, ^ 2- \ ^) (13) 

i=i i=i 

In the above formulas, the first equality shows the summation of the probabilities that index i is 
transmitted over the feedback channel, index j is received at the transmitter, and a pairwise error 
occurs. The second equality is the application of the Bayes' rule. Finally, the third one can be 
written noting that any two modulation matrices Zj and Z'j are in the constellation set j, chosen 
given index j at the transmitter. Hence, the event Zj ^ Zj \ i is conditioned on j. 

According to the average PEP expression in (fT3l) . the minimum PEP precoder set {Pi, ■ ■ ■ , Vk} 
can be obtained by solving the following optimization problem: 

min Ef=iEiiPf(jK)p(^) %,v. { {i+vch^v.pjny'' } 

s.t. Vj, VjV] y and Ti{VjV]) = N 

Here, VjVj >- shows that the Hermitian matrix VjVj is positive semi-definite. We have also 
omitted the constant term {1 + rjc) ~ from the objective function. Equation (fT4l) implies that 
deriving the set of precoder matrices depends on the feedback indices at the input/output of the 
quantizer, the set of Voronoi regions V, and the conditional index inversion probability of the 
feedback channel, Pi{j\i). 
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III. Quantizer and Precoder Codebook Design 

COVQs have been proposed to minimize the quantization average distortion when the quantiza- 
tion indices go through erroneous channels [35]- [38]. We use this idea to design a CSI quantizer 
for our system. Our COVQ design problem is to find the pair of Voronoi regions and precoder 
matrices (V, V) that solves the optimization problem expressed in (fT4l) . 

First, note that given a fixed set of Voronoi regions V, the joint optimization in ((14)) can be 
decoupled into a series of individual optimizations for each Vj as 

min E£iPf(jK)p(0 %ev. { (l + r/c ^V.PJn)"'^ | 
s.t. VjV] y and Tr(P,P]) = A^ 

This property simplifies the design procedure significantly. In ([T?]) . both the objective function and 
the constraints are convex functions of Vj. Furthermore, the objective function is differentiable 
throughout the whole domain of Vj . Therefore, we can use a steepest descent type of algorithm, 
such as the gradient algorithm to solve this problem [46]. We successively find the Voronoi regions 
associated to each index i from (O and the precoder matrices from separate implementations of 
([T?]) for every j. This algorithm proceeds as follows: 

A. Training-Based Gradient Algorithm: 

1) Generate a sequence of training samples of h by normalizing a sequence of A^-dimensional 
complex Gaussian random vectors. Then assume an initial set of positive semi-definite precoder 
matrices V = {Vi,--- , Vr} with powers \\Vj\\\ = iV, V j. 

2) Assign index i to each training vector h if 

i = arg min NJ Pf(jk) i^ + Vc^ '^rj^j ^^^^ 

The set of training vectors with assigned index i statistically represent Vj. Note that the above 
formula will be used later in the encoder of the quantizer after the codebook V is designed. 

3) To optimize the precoder matrix from ([T?]) . first define the following objective function for 
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each index j: 

J(^i) = Epf(jK)p(0 %ev. I (i + ^c h^v.vjhy'' I (17) 

The gradient of J('Pj) can be expressed as 

VJ(P,) = -2iVr/, J]pf(jK)p(^)E^,i,^ (^l + ^^^V,p]^)" ' hh^V, \ (18) 

In numerical implementation of (fTSi) . note that the random vector /i is uniformly distributed on 
the region Vj and its pdf is proportional to the volume of Vi, or the marginal probability p(z). For 
any function \1/, to find p(i) E^^y {^(/i)}, it is sufficient to add the values of '^{h) throughout 
the partition of the training space that is represented by index i (approximation of Vj) and divide 
the result by the size of the training sequence. 

4) To proceed with the gradient algorithm, update each matrix in the precoder codebook V 
using the following relation: 

n + 



V,{t + l) 



V,{t) - a{t) VJ(P,(t)) ^^ 



(19) 



In the above formula, Vj{t) denotes the value of the precoder matrix Vj in the t-th iteration of 
the algorithm. The step size a{t) can be set to any decreasing function of t, but from [46] we 
use (l+m)/(l+t) with an arbitrary positive real number m. The operation [ ■ ]^ shows that the 
matrix in the brackets is projected onto the space of positive semi-definite matrices with power A^. 
We use the Euclidean distance as the projection measure. With this measure, the above projection 
is simply removing all the negative eigenvalues of the matrix inside the brackets and normalizing 
the matrix, so that its power is N [47]. 

5) Use the output of the above algorithm and improve it by successively implementing steps (2) 
and (3-4) until convergence. Since the resulting sequence of objective function values from (fl?)) is 
decreasing and PEP is bounded bellow, the convergence of this algorithm is guaranteed. However, 
we cannot claim that by the above alternation we converge to a globally optimal solution. 

After storing the precoder matrix codebook V at the transmitter/receiver ends of the system, the 
encoder of the quantizer at the receiver side operates similar to Equation ([T6l) . At the transmitter 
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side, the task of the quantizer decoder is to choose Vj from the precoder codebook. 

In general, we can implement different index mapping methods in our scheme. For incorporating 
a mapping other than identity, we map the output index of the quantizer i to i' at the receiver side 
and pass it through the feedback channel. Then at the transmitter side, we use inverse mapping 
(demapping) to obtain the quantization index j from the output of the feedback channel, i.e., j' . 
Note that to drive the system with different mappings, one should replace Pf(j|i) with Pf(j'|i') in 
the design of the precoder codebook and in the implementation of the quantizer's encoder. 

B. System Design Without Knowing pf A Priori 

The design techniques established in the previous section can be extended to the case that the 
knowledge of pj is not accurate. Suppose that we know a coarse range of feedback bit error 
probability pf at the transmitter/receiver sides. The uncertainty about the feedback channel model 
can be taken into account assuming that fa ^ Pf ^ ft, where fa,fb G [0,0.5]. This type of 
channel modeling mismatch results from the uncertainty about the feedback channel conditions in 
a wireless environment and the uncertainty about the amount of resources that different receivers 
in the network may spend for protecting feedback signals. 

We develop our design strategy based on assuming a fixed pd E [fa, ft], which is called the 
design cross-over probability parameter. Inspired from [39], it is known that by designing COVQs 
based on worst-case assumptions, the system enjoys desirable robustness against channel modeling 
mismatch. Therefore, we find the set of precoder matrices and Voronoi regions that minimize the 
worst possible pairwise error probability. Similar to (fT4l) . the precoder design criterion can be 
expressed as: 

min max Ej^i E^^i Pd(j|OpW %ev, | (l + ^c ^^^i^]^) [ 

p, G [fajb] s.t. V J, P,P] y and Tr(P,PJ) = N 

where Pd(jK) is the conditional index inversion probability between i and j assuming that the 

BSC cross-over probability is p^. It is straightforward to show that the objective function in (|20l ) 

is an increasing function of the design cross-over probability parameter p^. We can show this 
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property by plotting Pc from (fT3] ) versus the SNR. Therefore, the min max value of pd coincides 
with the maximum point of the cross-over probability range, which is pd = maxpj = f^. By 
this choice, the rest of the COVQ design procedure from Section IIII-AI can be applied. 

IV. Numerical Results 

A. Precoder Characteristics 

In this section, we study the precoder structure of a 4-antenna system (M = 4) with 16 feedback 
indices {K = 16) or 4 feedback bits (log2 K = 4). The precoder codebook is designed for 
different cross-over probabilities pf. The solutions are obtained by solving optimization problem 
(O, using the training based gradient algorithm explained in Section IIII-AI Fig. [21 depicts the 
eigenvalues of the first member of the precoder codebook, Vi, i.e., [5i,- ■ ■ ,64]. Other members 
of the codebook also have similar properties. Here, ViVl = Ui Diag [5^, • • • , 5f\ U\, where Ui is 
unitary and Diag[-] denotes a diagonal matrix. Note that the transmitter's power constraint requires 
that YliK=i ^K = 4. When the feedback link is error-free or pf = 0, all the transmission power is 
assigned to the first eigenmode, denoted by the largest eigenvalue, 61. In other words, to achieve 
minimum PEP, transmission power is projected onto the direction of the major eigenvector of 
the precoder matrix, i.e. a column of Ui. This result coincides with the optimality of directional 
beamforming. With an error- free feedback link, we can design a directional beamforming system, 
where the transmit signals are the product of scalar symbols and unit- norm beamforming vectors. 
Hence, our precoding system degenerates to a beamforming system, where the first eigenvector 
of the precoder matrix is equivalent to the beamforming direction. As the cross-over probability 
of the feedback link increases, the eigenvalues of the minimal PEP precoder matrix converge to 
equal values. A better error performance can be obtained by spreading the power among different 
directions in space. In this case, our precoded matrix constellations also converge to matrices with 
equal eigenvalue squares, similar to open-loop orthogonal space-time block codes. 

B. BER Computations 

Fig.[3]shows the Bit Error Rate (BER) of our combining scheme over a 4x 1 MISO channel using 
Monte Carlo simulations. The feedback link carries 4 bits with different cross-over probabilities. 
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This situation resembles a high-rate feedback assumption K > M. The transmitter uses BPSK 
constellation symbols and the transmission rate is 1 bit/sec/Hz. In each data frame we transmit 130 
data symbols. In this figure, we demonstrate the performance of an open-loop STBC and that of the 
combining scheme without feedback errors. By increasing the cross-over probability pf from 0.00 
to 0.50, the BER of the system ranges between a directional beamforming system and an open- 
loop system. All the curves show full-diversity order within the range of SNRs that we consider 
in our simulations. This figure shows that even with feedback errors, our combining scheme 
preserves full-diversity order and it also provides additional array gain compared to an open-loop 
STBC. In this figure, we also plot the BER of Grassmannian beamforming from [15] [25] and 
show that it is optimal in terms of minimizing the PEP when feedback is error-free. According 
to our numerical experiments, different index mapping/demapping schemes at the input/output 
of the feedback channel result in similar performance results if the cross-over probability of the 
feedback link is known and the precoder codebook is optimized for each specific mapping. 

The BER performance of our scheme with low-rate and also erroneous feedback is illustrated in 
Fig. m In this experiment, we consider a 6 x 1 MISO channel with 2 feedback bits or 4 feedback 
regions. Directional beamforming in this case cannot achieve full-diversity order. Our combining 
system also shows the same property if we use 6-dimensional precoder matrices. In order to 
obtain full-diversity order, PSTBCs in the literature for instance the unitary PSTBCs in [25] use 
6x2 precoder matrices from unitary constellations in [48], applied to 2 x T orthogonal design 
matrices [42]. We use generalized POD constellation matrices similar to the one introduced in 
([3]), based on a 6 x 8 orthogonal code structure. This code is generated by removing two rows of 
an 8 X 8 orthogonal design, leaving the first two rows without precoding, and multiplying 4x4 
precoder matrices inside the code structures at the lower 4 rows of the matrix. The simulation 
curves can be extended using other values of A^, which we skip for the sake of clarity in the figure. 
Our conclusion from this study is that only when N < K, full diversity order can be obtained. 
Moreover, when K < M, minimum BER performance is attributed to the POD structure with 
ii"-dimensional precoders. Over an error- free feedback channel, our combining scheme with 4- 
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dimensional precoders outperforms PSTBCs with 6x2 unitary precoders. Fig. |4] also shows the 
performance of the system with different feedback cross-over probabilities. With a non-zero cross- 
over probability in the feedback link, PSTBCs with unitary precoders fail to obtain full-diversity 
order. Our combining scheme with 4-dimensional precoder matrices, designed based on knowing 
Pf, however preserves full-diversity order, even if feedback is low-rate and erroneous. 

In Fig. [5l we show the performance of our combining schemes over a 4-antenna MISO channel 
with complex constellations. There is no M x T rate-one complex orthogonal design for M > 2 
[42]. In this case, instead of using orthogonal design structures with less rates, we can use a 
quasi-orthogonal space-time block code (QOSTBC) [43] [44], to obtain a generalized partly 
quasi-orthogonal design (PQOD) structure. The quasi-orthogonal inner code of this example 
uses QPSK modulation symbols with j rotations. The transmission rate is 2 bits/sec/Hz. With 
K = A and K = 2 feedback regions, we use generalized PQODs with A^ = 4 and N = 2- 
dimensional precoders, respectively. Again, with feedback errors, these constellations can provide 
full-diversity order and additional array gain compared to open-loop QOSTBCs. In this figure, 
we also demonstrate the performance of this system employing Grassmannian beamforming, 
when feedback is erroneous. This situation resembles a noiseless design mismatch and results 
in not achieving full-diversity order. Here, the number of feedback regions is equal to the number 
of transmit antennas and Grassmannian beamforming is similar to antenna selection (AS). The 
performance of our combining scheme designed assuming noiseless feedback (VQ-design) is 
also very close to that of Grassmannian beamforming. In this case, there are N independent 
beamforming directions in space with similar cross (Chordal) distances [15], [25]. Upon receiving 
an erroneous feedback index, it does not matter which direction we pick. Hence, employing 
different index mapping schemes in this case does not change the performance of the system, 
even when the feedback channel is erroneous. 

C. Effect of Mismatch 

In the last experiment, we investigate the effects of channel modeling mismatch on the perfor- 
mance of our system. For the numerical results shown in Fig. [6l we again consider a 4 x 1 MISO 
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channel with 16 feedback regions and BPSK symbol constellations. We assume certain feedback 
cross-over probabilities for designing the COVQ or VQ-based combining system and use different 
ones for Monte Carlo simulations. The following cases are studied: 

First, the precoder codebook is designed assuming feedback is noiseless. We call this design 
a VQ-based one since COVQ becomes a VQ when p^ = 0.00. Also the combining scheme 
degenerates to directional beamforming. The resulting system is examined using a feedback link, 
where the actual cross-over probability is p/ = 0.04. The performance results are very close 
to the performance of Grassmannian beamforming when operating over an erroneous feedback 
channel with pf = 0.04. This type of mismatch, which results in severe performance degradation 
is referred to as VQ-based design mismatch. In this case, since the number of feedback regions 
is more than the number of transmit antennas, different beamforming directions have different 
cross (Chordal) distances [15], [25]. Hence, we can improve the performance of the system by 
modifying the feedback index mapping without any additional cost. To show this property, we first 
define the average Chordal distance of the beamforming directions at the transmitter, conditioned 
on the feedback indices at the receiver, when the feedback cross-over probability is given. Then we 
minimize this average distance over the space of possible index mapping solutions using simulated 
annealing (SA) [37]. The details of this algorithm are skipped for brevity. The resulting VQ-based 
system outperforms the same system with identity mapping. However, the gain of optimizing the 
index mapping is not significant compared to the gain of COVQ-based design since it does not 
change the diversity of the system. Second, we assume the design cross-over probability pd = 0.04 
and examine the system performance over an error-free feedback link, where pj = 0.00. This 
situation is called COVQ-based design mismatch. This type of mismatch does not degrade the 
performance of the system severely. The performance results in this case are worse than those 
of a system with ideal design of precoders {pd = 0.00) and error-free simulations (p/ = 0.00). 
However, full-diversity order is still preserved. 

Third, we consider a feedback channel where the cross-over probability changes uniformly in 
the range 0.00 < p/ < 0.04. In this scenario, as we explained in Section ITlI-Bi we use a worst-case 
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design strategy to obtain robustness against feedback channel modeling mismatch. The precoder 
codebook is designed assuming that the cross-over probability is p^ = 0.04. By this choice, full- 
diversity order is preserved and also the array-gain of the system is superior to an open-loop 
system. Fourth, we consider the same feedback channel model with the same range of cross-over 
probabilities, but this time we assume an alternative design parameter, where p^ is the average 
value of feedback errors, i.e., pd = (fa + /fe)/2 = 0.02. By this assumption, the transmission 
scheme does not achieve full-diversity order. The latter observation confirms that the worst case 
design assumption pd = max^j-^ j^] pf is a better candidate compared to the average error design. 

V. Summary and Conclusions 

In this paper, we studied high-rate and low-rate feedback (closed-loop) communication systems 
with possible feedback errors. We showed numerically that our combining schemes obtain full- 
diversity order with additional array gain compared to open-loop systems. The design criterion 
of our scheme was derived based on a pairwise error probability measure and our system was 
optimized using a training-based algorithm. As the feedback cross-over probability approaches 
zero, our scheme converges to a directional beamforming system. On the other hand, as the 
feedback error increases, our scheme converges to a no-CSIT or open-loop system. We presented 
combining strategies for both high-rate and very low-rate feedback scenarios, which is simply 
extendable to the cases that the feedback link is noisy. With very low-rate feedback, the dimension 
of the precoder matrix in our scheme reduces and our system converges to an open-loop STBC 
structure. 

The design strategies and solution algorithms presented in this paper are robust against feedback 
channel modeling mismatch. Using a worst-case design approach, we designed the quantizer 
structure and derived the precoder matrices used in the code structures, even if the exact value 
of the feedback channel error is unknown to the transmitter/receiver sides a priori. It was shown 
numerically that even in this case, our transmission scheme preserves full-diversity order and 
provides additional array gain compared to open-loop systems. 



19 

References 

[I] IEEE Std 802. 16e^^-2004/Corl-2005, IEEE standard for local and metropolitan area networks, Part 16: air interface for 
fixed broadband wireless access systems, Ammendment 2: physical and medium access control layers for combined fixed and 
mobile operation in licensed bands, and corrigendum 1, February 2006. 

[2] A. Narula, M.L. Lopez, M.D. Trott, and G.W. Wornell, "Efficient use of side information in multiple- antenna data transmission 

over fading sliannels," IEEE Journal on Selected Areas in Communications, vol. 16, October 1998, pp. 1423-1436. 
[3] S.A. Jafar and A.J. Goldsmitfi, 'Transmit optimization and optimality of beamforming for multiple antenna systems with 

imperfect feedback," IEEE Trans, on Wireless Communications, vol. 3, July 2004, pp. 1165-1175. 
[4] M. Skoglund and G. Jongren, "On the capacity of multiple-antenna communication link with channel side information," IEEE 

Jour on Sel. Areas in Communications, vol. 21, Apr 2003, pp. 395-405. 
[5] Y. Xie, C. N. Georgiades, and A. Arapostathis, "Minimum outage probability transmission with imperfect feedback for MISO 

fading channels,", IEEE Trans, on Wireless Comm, vol. 4, May 2005, pp. 1084-1091. 
[6] S. Bhashyam, A. Sabharwal, and B. Aazhang, "Feedback gain in multiple antenna sysmtes," IEEE Trans, on Communications, 

vol. 50, May 2002, pp. 785-798. 
[7] G. Jongren, M. Skoglund, and B. Otterson, "Combining beamforming and orthogonal space-time block coding," IEEE Trans. 

on Information Theory, vol. 48, March 2002, pp. 611-627. 
[8] S. Zhou and G.B. Giannakis, "Optimal transmitter eigen-beamforming and space-time block coding based on channel mean 

feedback", IEEE Trans, on Signal Processing , vol. 50, October 2002, pp. 2599-2613. 
[9] S. Zhou and G.B. Giannakis, "Optimal transmitter eigen-beamforming and space-time block coding based on channel 

correlations," IEEE Trans, on Information Theory, vol. 49, July 2003, pp. 1673-1690. 
[10] H. Sampath and A. Paulraj, "Linear precoding for space-time coded systems with known fading correlations," IEEE 

Communications Letters, vol. 6, June 2002, pp. 239-241. 

[II] L. Liu and H. Jafarkhani, "Application of quasi-orthogonal space-time block codes in beamforming," IEEE Trans, on Signal 
Processing, vol. 53, January 2005, pp. 54-63. 

[12] V. Lau, Y. Liu, and T.-A. Chen, "On the design of MIMO block-fading channels with feedback-link capacity constraint," 

IEEE Trans. Communications, vol. 52, Jan. 2004, pp. 62-70. 
[13] Y. Ko and C. Tepedelenlioglu, "Orthogonal space-time block coded rate-adaptive modulation with outdated feedback," IEEE 

Trans, on Wireless Communications, vol. 5, February 2006, pp. 290-295. 
[14] K.K. Mukkavilli, A. Sabharwal, E. Erkip, and B. Azhang, "On beam forming with finite rate feedback in multiple-antenna 

systems,", IEEE Trans, on Information Theory, vol. 49, October 2003, pp. 2562-2579. 
[15] D.J. Love, R.W. Heath, and T. Strohmer, "Grassmannian beamforming for multiple-input multiple-output wireless systems," 

IEEE Trans, on Information Theory, vol. 49, October 2003, pp. 2735-2747. 
[16] D.J. Love and R.W Heath, "Limited feedback unitary precoding for spatial multiplexing systems," IEEE Trans, on Information 

Theory, vol. 51, August 2005, pp. 2967-2976. 
[17] J.C. Roh and B.D. Rao, "Transmit beamforming in multiple antenna systems with finite rate feedback: a VQ-based approach," 

IEEE Trans, on Information Theory, vol. 52, March 2006, pp. 1101-1112. 
[18] J. Zheng, E. Duni, and B. D. Rao, "Analysis of multiple-antenna systems with finite-rate feedback using high resolution 

quantization theory," IEEE Trans, on Signal Processing, vol. 55, pp. 1461-1476, April 2007. 
[19] S. Zhou, Z. Wang, and G.B. Giannakis, "Quantifying the power loss when transmit beamforming relies on finite rate 

feedback," IEEE Trans, on Wirel. Communications, vol. 4, July 2005, pp. 1948-1957. 
[20] Y Li, N.B. Mehta, A.F. Molisch, and J. Zhang, "Optimal signalling and selection verofocation for single transmit antenna 

selection," IEEE Trans, on Communications, vol. 55, April 2007, pp. 778-789. 
[21] A. Hottinen, O. Trikkonen, and R. Wichman, "Multi-antenna transciever techniques for 3G and beyond," John Willey, 2003. 
[22] H. S. Wang and N. Moayeri, "Finite-state Markov channel- a useful model for radio communication channels," IEEE Trans. 

on Vehicular Technology, vol. 44, February 1995, pp. 163-171. 
[23] S. Ekbatani, F Etemadi, and H. Jafarkhani, "Transmission over slowly fading channels using unreliable quantized feedback," 

IEEE Data Compression Conf, March 2007. 



20 

[24] K. Huang, B. Mondal, R.W. Heath, and J.G. Andrews, "Effect of feedback delay on multi-antenna limited feedback for 

temporally-correclated channels," IEEE Global Telecom. Conference, November 2006. 
[25] D.J. Love and R.W. Heath, "Limited feedback unitary precoding for orthogonal space-time block codes," IEEE Tram, on 

Signal Processing, vol. 53, January 2005, pp. 64-73. 
[26] Y. Yu, S. Keroueden, and J. Yuan, "Closed-loop extended orthogonal space-time block codes for three and four transmit 

antennas," IEEE Signal Processing Letters, vol. 13, May 2006, pp. 273-276. 
[27] D.A. Gore and A.J. Paulraj, "MIMO antenna subset selection with space-time coding," IEEE Trans, on Signal Processing, 

vol. 50, October 2002, pp. 2580-2588. 
[28] A. Zhou, J. Wu, Z. Wang, and M. Doroslovacki, "Quantifying the performance gain of direction feedback in a MISO system," 

IEEE International Science and Systems Conf., March 2006, pp. 573-578. 
[29] E.G. Larsson, G. Ganesan, R Stoica, and W.-H. Wong, "On the performance of orthogonal space-time block coding with 

quantized feedback," IEEE Communications Letters, vol. 6, November 2002, pp. 487-489. 
[30] J. Akhtar and D. Gesbert, "Extending orthogonal block codes with partial feedback," IEEE Trans, on Wireless Communica- 
tions, vol. 3, November 2004, pp. 1959-1962. 
[31] J.K. Milleth, K. Giridhar, and D. Jalihal, "Pre-processed space-time trellis codes with one-bit feedback," IEEE Communica- 
tions Letters, vol. 9, Aug 2005, pp. 703-705. 
[32] L. Liu and H. Jafarkhani, "Space time trellis codes based on channel phase feedback, " IEEE Trans, on Communications, 

vol. 54, December 2006, pp. 2186-2198. 
[33] G. Jongren and M. Skoglund, "Quantized feedback information in orthogonal space-time block coding," IEEE Trans, on 

Information Theory, vol. 50, October 2004, pp. 2473-2482. 
[34] A. Papoulis and S.U. Pillai, "Probability, random variables and stochastic processes," McGraw Hill, 4th Edition, 2002. 
[35] H. Kumazawa, M. Kasahara, and T. Namekawa, "A construction of vector quantizers for noisy channels," Electronics Eng. 

in Japan, vol. 67 -B, January 1984. 
[36] N. Farvardin and V. Vaishampayan, "Optimal quantizer design for noisy channels: an approach to combined source-channel 

coding," IEEE Trans, on Information Theory, vol. 33, November 1987, pp. 827-838. 
[37] N. Farvardin, "A study of vector quantization for noisy channels," IEEE Trans, on Information Theory, vol. 36, July 1990, 

pp. 799-809. 
[38] N. Farvardin and V. Vaishampayan, "On the performance and complexity of channel-optimized vector quantization," IEEE 

Trans, on Information Theory, vol. 37, January 1991, pp. 155-160. 
[39] H. Jafarkhani and N. Farvardin, "Design of channel-optimized vector quantizers in the presense of channel mismatch," IEEE 

Trans, on Communications, vol. 48, January 2000, pp. 118-124. 
[40] S. Ekbatani and H. Jafarkhani, "Combining beamforming and space-Time coding using quantized feedback," In revision for 

IEEE Trans, on Wireless Communications. http://www.ece.uci.edu/~sekbatan/Twir_SEHJ_8.pdf 
[41] V. Tarokh, H. Jafarkhani, and A.R. Calderbank, "Space-time block codes from orthogonal designs," IEEE Trans, on 

Information Theory, vol. 45, July 1999, pp. 1456-1467. 
[42] H. Jafarkhani, "Space-Time Coding: Theory and Practice," Cambridge University Press, 2005. 
[43] H. Jafarkhani, "A Quasi-Orthogonal Space-Time Block Code," IEEE Trans, on Communications, vol. 49, no. 1, Jan 2001, 

pp. 1-4. 
[44] W. Su and X-G Xia, "Signal constellations for quasi-orthogonal space-time block codes with full-diversity," IEEE Trans, on 

Information Theory, vol. 50, October 2004, pp. 2331-2347. 
[45] S. Ekbatani and H. Jafarkhani, "Design of multi-antenna coded modulators using noisy quantized channel state information," 

IEEE Global Telecom. Conf, November 2006. 
[46] D.R Palomar and M. Chaing, "A tutorial on decomposition methods for network utility maximization," IEEE Jour on Sel. 

Areas in Communications, vol. 24, August 2006, pp. I439-145I. 
[47] S. Boyd and L. Vandenberghe, Convex Optimization, Cambridge University Press, 2004. 
[48] B. M. Hochwald, T. L. Marzetta, T J. Richardson, W. Sweldens, and R. Urbanke, "Systematic design of unitary space-time 

codnstellations," IEEE Trans, on Information Theory, vol. 6, September 2000, pp. 1962-1973. 



21 



Transmitter 



POD'S 
inner STBC 


— 


POD'S 
precoder 












Input 
Bits 


COVQ 

Decoder 



Fading | 

■ 

Channel ■ 



Fig. 1. System block diagram. 



Pf(7l0 



Feedback Link 



Receiver 



Channel Estimation 
and Synchronization 



COVQ 

Encoder 



ML 

Decoder 



Output 
Bits 





- q 








2 


-O- Pj = 0.00 _ 




\ 






^p, = 0.04 


1.8 


\ 
\ 
\ 






-^P, = 0.10 
^-P, = 0.40 


1.6 


- 


1.4 


- ^^v\ 






- 


1.2 


\ 

\ 


^^^^^ 




- 






- ^ J^^XN|^ 




~ ~ -V 




^X:^^ 




0.8 


- 


x. 


^"^"~^^ 


- 


0.6 


- 


\ 
\ 
\ 


^\x^""^^^\ 


- 


0.4 


- 


\ 
\ 


N^^^ 


"^-^ - 


0.2 


- 


\ 
\ 




^^ 





1 


1 




o 

1 



Fig. 2. Eigenvalues of "Pi, [<5i, ■ ■ ■ , ^4] for an Af = 4 multi-antenna system with K — 16 feedback regions and different 
cross-over probabilities. 



22 



LU 

m 




-^- 4-ant, open-loop OSTBC 
„ POD, 4-bit COVQ, p = p^= 0.40 

^ POD, 4-bit COVO, pj= p^= 0.10 

A. POD, 4-bit COVO, pj= p^= 0.04 

_Q_ POD, 4-bit VQ, p = p^= 0.00 

Grassmann. VQ, Pj= 0.00 



10 
SNR(dB) 



11 



12 



13 



14 



15 
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